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Abstract 

The electromagnetic form factors of the trinucleons 3 H and 3 He are calcu- 
lated with wave functions obtained with the Argonne v\s two-nucleon and 
Urbana IX three-nucleon interactions. Full account is taken of the two-body 
currents required by current conservation with the v\s interaction as well as 
those associated with NA transition currents and the currents of A resonance 
components in the wave functions. Explicit three-nucleon current operators 
associated with the two-pion exchange three-nucleon interaction arising from 
irreducible S-wave pion-nucleon scattering is constructed and shown to have 
very little effect on the calculated magnetic form factors. The calculated mag- 
netic form factor of 3 H, and charge form factors of both 3 H and 3 He are in 
satisfactory agreement with the experimental data. However, the position of 
the zero in the magnetic form factor of 3 He is slightly underpredicted. 

21.45.+v,25.30.Bf,27.10.+h 
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I. INTRODUCTION 



The electromagnetic form factors of the few-body nuclei, along with the deuteron struc- 
ture functions and cross section for threshold electrodisintegration at backward angles, are 
the observables of choice for testing the quality of models for the nuclear interaction and 
the associated current operator, including its exchange current components. Such testing 
has become possible by the development of practical computational methods for numerical 
calculation of the wave functions of the few nucleon systems, which correspond to realistic 
phenomenological interaction models [1]. Employing such wave functions along with the 
the two-nucleon exchange current operators, which are required by the continuity equation 
and/or consistency with the interaction models, e.g., by Poincare invariance, it has become 
possible to predict the experimental electron scattering observables of the few nucleon sys- 
tems up to momentum transfers of about 2 GeV/c in an at least qualitatively satisfactory 
way. 

Among the remaining open issues are the need for quantitative understanding of the 
form factors of the trinucleons in the region around and above their first zeros. While the 
long standing, unsettled issue of the behavior of the tensor polarization of the deuteron for 
momentum transfers above 3 fm -1 appears close to settlement by high quality experimen- 
tal work at the Thomas Jefferson National Accelerator Facility, there remains a need for 
potential model development in the case of the trinucleon systems. This is partly because 
of the remaining problem in quantitative understanding of the trinucleon form factors at 
high momentum transfer and partly because of the resilient open issues concerning the form 
of the three-nucleon interaction, which appears to be required for the understanding of the 
binding energies of the light-nuclei (A <7) ground states @. 

Here this question is investigated in several different ways. First, a numerically extensive 
calculation of the electromagnetic form factors of the trinucleons is presented with high pre- 
cision variational wave functions, which correspond to the Argonne v is two-nucleon || and 
the Urbana IX three-nucleon [|J interactions. In this calculation the two-nucleon exchange 
current operators are constructed by the same method as used in the earlier calculations 
in Refs. that employed the Argonne v u interaction J/J. Second, the irreducible three- 
nucleon exchange current operator, which corresponds to the best understood part of the 
two-pion exchange three-nucleon interaction associated with 5-wave pion nucleon scattering 
on the intermediate nucleon is constructed, and its matrix elements for the trinucleon bound 
states are shown to be very small for momentum transfer values below 1 GeV/c. Third, a 
systematic treatment of A-isobar configurations in the trinucleon ground states is made, and 
their effect on the trinucleon form factors are calculated with inclusion of all the associated 
and required exchange current operators. 

The calculated magnetic form factor of 3 H is in fairly good agreement with the present 
experimental data once the exchange current contributions are included. That of 3 He agrees 
less well with the corresponding data at high values of momentum transfer, which is a 
consequence of the underpredicted position of the first zero in the form factor (it falls at 
3.75 fm -1 , which is below the experimental range 4.2-4.4 fm -1 ). This result is independent 
of the presence or absence of the A-isobar configurations in the wave function model. The 
problem is likely to have its origin in a somewhat too weak overall strength of the model 
for the isovector part of the exchange current operator at large momentum transfer. As 
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this is constructed so as to be consistent with the vi$ two-nucleon interaction, the ultimate 
origin of the problem with the high momentum behavior of the calculated isovector magnetic 
form factor of the trinucleons may reside with the potential model if not with some purely 
transverse exchange current mechanism that has not been considered. The effect of the 
irreducible three-nucleon exchange current operator on this form factor is very small. The 
finding that the A-isobar configurations have a very small effect on the trinucleon form 
factors conforms to earlier results obtained with other potential models [^]. 

The calculated charge form factor of 3 He agrees very well with the experimental values 
over the measured range of momentum transfer, with an exception for its highest end. In 
the case of 3 H quantitative agreement with the experimental form factor is achieved only up 
to the position of the secondary maximum, above which region the calculated values are too 
large by factors 2-3. The exchange current contributions are essential for agreement with 
the experimental charge form factors. 

This paper is divided into four main sections. Section II contains a description of the 
calculation of the electromagnetic form factors of the trinucleons using a purely nucleonic 
variational wave function constructed for the Vi$ model augmented by the Urbana IX three- 
nucleon interaction. The details of the hyperspherical variational model wave function are 
given in Subsection II-A. Subsection II-B contains the description of the model for the elec- 
tromagnetic current operator including the exchange current operators. In Subsection II-C 
the irreducible three-nucleon exchange current operator, which corresponds to the main 
nonresonant two-pion exchange three-nucleon interaction, is derived. Finally Subsection 
II-D contains the form factor results obtained with this restricted model. In Section III 
the description of the extended model wave function and current operators that include the 
A-isobar configurations in the wave function is given. The model for the NA transition 
potential is described in Subsection III-A and the corresponding current operators are de- 
scribed in Subsection III-B. The calculation of the form factors in the extended model is 
outlined in Subsections III-C and III-D. Finally Section IV contains a concluding discussion. 

II. TRINUCLEON FORM FACTORS WITH NUCLEONIC WAVE FUNCTIONS 

In this section the calculation of the elastic form factors of the A=3 nuclei with wave 
functions for a realistic Hamiltonian model formed of the Argonne V\g (AV18) two-nucleon || 
and Urbana IX (UIX) three-nucleon Q interactions is described. The calculation employs 
charge and current operators that besides the standard single nucleon components also 
contain two-nucleon components, the leading terms of which are constructed consistently 
with the AV18 model. The three-nucleon exchange current operator, which corresponds to 
the two-pion exchange three-nucleon interaction associated with isospin odd S-w&ve pion 
rescattering is also derived and is shown to have only a minor effect on the trinucleon form 
factors. 

A. The AV18/UIX Hamiltonian and trinucleon wave functions 

The AV18 model [[J is a recent high-quality nucleon-nucleon interaction containing ex- 
plicit charge-symmetry-breaking (CSB) and charge-independence-breaking (CIB) terms, as 
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well as a complete treatment of the electromagnetic interaction up to order a 2 , a being the 
fine structure constant. It is constructed to fit the Nijmegen pp and np scattering database, 
low-energy nn scattering parameters, and the deuteron binding energy with a x 2 P er datum 
close to 1. 

The UIX three-nucleon interaction [|J consists of a long-range term due to excitation of 
an intermediate A-isobar via two-pion exchange and a short-range repulsive phenomenologi- 
cal term, which simulates the dispersive effects which arise upon integrating out A-degrees of 
freedom. The strength of this repulsive term is determined by fitting the triton binding en- 
ergy in "exact" Green's Function Monte Carlo (GFMC) calculations [[J and the equilibrium 
density of nuclear matter in variational calculations based on operator-chain summation 
techniques ||. 

Recent GFMC calculations based on the AV18/UIX Hamiltonian model have been shown 
to provide a good description of the low-energy spectra and charge radii of nuclei with 
A <7 0. In particular, the calculated binding energies of 3 H and 3 He are within a few keV 
of the experimental values (see Table |). 



In the present work we use trinucleon wave functions obtained by Kievsky et al. ||TQ , |lT| . |T2 
with the pair-correlated hyperspherical harmonics (PHH) method. Although variational this 
method has been refined in the last few years so much that it yields results with an accuracy 
comparable to that achieved in recent Faddeev and GFMC calculations, as may be seen 
in Table |. The PHH method is briefly reviewed here, for completeness; however, a more 
thorough discussion of it as well as its extensions to describe both the A=3 low energy 



continuum and A=A ground state can be found in Refs. fllT[|l4|13 . 

The wave function \l/ of a three-nucleon system with total angular momentum JJ Z and 
total isospin TT Z can be decomposed as 

¥ = X>(Xi,y<) , (2.1) 

where the amplitude ?/>(xj, y^) is a function of the Jacobi coordinates X; = — and 
yi = {vj + v k — 2rj) /y/3, i,j,k being a cyclic permutation of 1,2,3. To ensure the overall 
antisymmetry of the amplitude T/>(xj, y*) is antisymmetric with respect to exchange of 



nucleons j and k. In the PHH method, it is expressed as [11,12] 



V>(xi, Yi) = fa(xi)^ a (xi, yi)y a {j, k; i) , (2.2) 



y a (j,k; 
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J A a 









TT Z 



(2.3) 



L a and A a , the spin 



where each channel a is specified by the orbital angular momenta £ a , 
(isospin) S£ (T^ k ) of pair jk and the total spin S a . Orbital and spin angular momenta are 
coupled, in the LS-scheme, to give total angular momenta JJ Z . The (channel-dependent) cor- 
relation functions f a (xi) are obtained from solutions of two-body Schrodinger-like equations 



in channel jplpSp |I0|J11|1 , and take into account the strong state-dependent correla- 
tions induced by the nucleon-nucleon interaction. They improve the behavior of the wave 
function at small interparticle distances. Were it not for their presence, the decomposition 
in Eq. (|2.2j ) would be identical to that in the Faddeev scheme J14 . 
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Next, the hyperspherical coordinates p and fa, defined as 

P = \/ x i+V? , cos (pi = Xi/p (2.4) 
are introduced, and the dependence of & a (xi,yi) on p and fa is made explicit by writing 



(2.5) 



ra=0 



n 



COS( 



sin i 



p 



cos 20j 



(2.6) 



where N^ a,La are normalization factors, P"^ are Jacobi polinomials and n is a non-negative 
integer, n = 0, ■ ■ • , M a , where M a is the selected number of basis functions in channel a. 
The complete wave function is then written as 



ij'fc cyclic Q! n=0 

The Rayleigh-Ritz variational principle, 

(5 U V\H-E\V) = 0, 



(2.7) 



(2.8) 



is used to determine the hyper-radial functions u"(p) in Eq. (|2.7|). Carrying out the variation 
5 u ty with respect to the functions u%(p), the following equation is easily derived: 



E (fa{xi)y a {j, k; i)Z"(fa)\H - E\*&) | n = 

r/'fc cyclic 



(2.9) 



where f2 denotes the angular variables 0j, Xj and y*- Performing the integration over Q 
and spin-isospin sums (as implicitly understood by the notation (• • •) |q) leads to a set of 
coupled second order differential equations for the u®{p), which is then solved by standard 
numerical techniques [JTOlJTl . 



The binding energy of the A=3 nuclei obtained with the PHH method from the 
AV18/UIX Hamiltonian are listed in Table | [EJ. Also listed in Table [I] are results cal- 



culated with converged r-space |15j and p-space JL61 Faddeev wave functions for an older 
model of the two-nucleon interaction, the Argonne (AV14) 0. The binding energies 
obtained with the various methods are in excellent agreement with each other, typically 
within 10 keV or less. 



B. Nuclear charge and current operators 

A fairly complete decription of the model for the nuclear electromagnetic current has 
been most recently given in Ref. ]1|]. Here we only review its general structure. The nuclear 
charge and current operators are expanded into a sum of one- and two-body terms: 
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p(q) = £Pi (q)+£/f(q), (2-10) 

j(q) = Ej? ) (q) + Ejg ) (q)> (2.11) 

where q is the momentum transfer. The one-body operators and p 1 ) are given by 



q-fT.xp^, (2.12) 



(2.13) 

up to terms proportional to 1/m 2 , m being the nucleon mass. Equation ( |2.12| ) includes the 
leading relativistic corrections to single-nucleon charge operator, namely the Darwin-Foldy 
and spin orbit terms. Here the G S J.J M (q 2 ^) are the electric/magnetic (E/M) isoscalar/isovector 
(S/V) form factors of the nucleon, taken as function of the four-momentum transfer 

ql = q 2 -^ 2 >0 , (2.14) 

where the energy transfer uj = \J<i 2 + m 2 ^ — for elastic scattering on a target of mass 
ttit initially at rest in the lab. These form factors are normalized as 

G|(0) = G V E (0) = 1 , 
G S M (0) = 0.880 n N , 

G V M {Q) = 4.706 /ijv, (2.15) 

Hn being the nuclear magneton, and their ^-dependence is constrained by analyzing 
electron-proton and electron-deuteron scattering data. While the proton electric and mag- 
netic form factors are experimentally fairly well known over a wide range of momentum 
transfers, there is significant uncertainty in the neutron form factors, particularly the electric 
one, which are obtained from model-dependent analyses of ed data. Until this uncertainty in 
the detailed behaviour of the electromagnetic form factors of the nucleon is narrowed, quan- 
titative predictions of electro-nuclear observables at high momentum transfers will remain 
rather tentative. We will re-examine this issue in Sec. II-D-1 below. 



p, (1) (q) 



MqJ 



1 + ql/Am 2 2 



G s E (ql) + G v E (ql)T s 



8m 2 



2GM) - G S E {ql) + \2G\Ml) ~ Gl{ql) 



Am 



G E (ql) + G V E (ql)r z ,] {p*,^} 



4m 



GUQi) + G V M (ti)T; 



q x (Tie 



1. The two-body current operator 

The two-body current operator can be separated into a model-independent (MI) term 
determined from the interaction (in the present case, the charge-independent part of the 
AVI 8 model) following a prescription originally proposed in Ref. |I7J], and a model-dependent 
(MD) one, associated with the p7T7 and uttj electromagnetic couplings. Explicit expressions 
for all these currents have been most recently given in Ref. ]18| . 
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The p7T7 and ooir^y MD currents are purely transverse and therefore unconstrained by 
the nucleon-nucleon interaction. The values of the transition form factors G piT1 (q 2 ^) and 
Gr W7r7 (g2) at the photon point are known to be G p7ry (0) = g pnJ = 0.56, Ref. [Tj|, and 
C W7r7 (0) = g W7r7 = 0.68, Ref. p0| , from the measured widths of the p — > 717 and u; — > 7r7 
decays, while their ^-dependence is modeled using vector-meson dominance. Monopole 
form factors at the pion and vector-meson strong interaction vertices are introduced to take 
into account the composite nature of nucleons and mesons. The cutoff parameters A„-, A p 
and A w in these form factors are not known. Here we use the values A n = 0.75 GeV and 
A p = A w = 1.25 GeV obtained from studies of the B-structure function of the deuteron |f2T|. 

The leading MI two-body currents, denoted as pseudoscalar (PS) or 7r-like and vector 
(V) or p-like, are the isovector ones associated with the isospin-dependent central, spin-spin, 
and tensor components of the interaction. Their derivation has been given in a number of 
references |I|j22|, and will not be repeated here. We only note that: i) the PS and V two- 
body currents have no free parameters and, by construction, satisfy the continuity equation 
with the given realistic interaction (here the charge-independent part of AV18 model); ii) the 
continuity equation requires the same form factor be used to describe the electromagnetic 
structure of the hadrons in the longitudinal part of the current operator and in the charge 
operator, while it places no restrictions on the electromagnetic form factors which may be 
used in the transverse parts of the current. Ignoring this ambiguity, the form factor G^(g^) 
is used in the PS and V currents operators, in line with the "minimal" requirements of 
current conservation. 

There are additional two-body currents associated with the momentum dependence of 
the interaction, but their construction is less straightforward. A procedure similar to that 
used to derive the PS and V currents has been generalized to the case of the currents 
from spin-orbit components of the interaction It consists, in essence, of attributing 

these to exchanges of er-like and a>-like mesons for the isospin-independent terms, and to 
p-like mesons for the isospin-dependent ones. The explicit form of the resulting currents, 
denoted as SO, can be found in Refs. fl8|j23| . The two-body currents from the quadratic 



momentum dependence of the interaction are obtained by minimal substitution — ► — 
\ Ge(Qu) + Gg(gJ)r 2i j A(rj), A(r^) being the vector potential, into the corresponding 
components. In the case of the AV18 model, the p 2 -dependence is via L 2 and (L • <T\ L • cr 2 + 
h.c.) terms, and the associated currents are denoted respectively as LL and S02 fP8|,|5|. 



We note that the SO, LL and S02 currents are fairly short-ranged, and have both 
isoscalar and isovector terms. Their contribution to isovector observables is found to be 
numerically much smaller than that due to the leading PS (7rTike) current. However, these 
currents give non-negligible corrections to isoscalar observables, such as the deuteron mag- 
netic moment and B-structure function [p^j . Finally it is worth enphasizing that, while the 
construction in Ref. is not unique, it has nevertheless been shown to provide, at low and 
moderate values of momentum transfer, satisfactory description of most observables where 
the isovector two-body currents play a large (if not dominant) role, such as the deuteron 
threshold electrodisintegration the neutron and proton radiative captures on protons 



and deuterons at low energies |18|j2^| , and the magnetic moments and form factors of the 
trinucleons (as shown below). 
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2. Two-body charge operators 



While the MI two-body currents are linked to the form of nucleon-nucleon interaction via 
the continuity equation, the most important two-body charge operators are model dependent 
and may be viewed as relativistic corrections. They fall into two classes. The first class 
includes those effective operators that represent non-nucleonic degrees of freedom, such 
as nucleon-antinucleon pairs or nucleon-resonances, and which arise when these degrees of 
freedom are eliminated from the state vector. To the second class belong those dynamical 
exchange charge effects that would appear even in a description explicitly including non- 
nucleonic excitations in the state vector, such as the pn^ and uiry transition couplings. The 
proper forms of the former operators depend on the method of eliminating the non-nucleonic 
degrees of freedom E3,EB|f2^| . There are nevertheless rather clear indications for the relevance 



of two-body charge operators from the failure of calculations based on the one-body operator 
in Eq. (2.12) in predicting the charge form factors of the three- and four-nucleon systems ||, 



and deuteron A-structure function and tensor polarization observable [23,28 



The two-body model used in the present work consists of the 7T-, p- and a;-meson ex- 
change charge operators, as well as of the pir'y and utfy charge transition couplings. The 
former are derived by considering the low-energy limit of the relativistic Born diagrams 
associated with the virtual meson photoproduction amplitude. The pir'y and um'y opera- 
tors are the leading corrections obtained in a non- relativistic reduction of the correspond- 
ing Feynman diagrams with transition couplings, for example (ir a (k) | j M (0) | p b (p, e)) = 
— G p7r7 (g^)/m p dabt^var P u k a e T , e being the polarization vector of the p-meson. Coupling 
constants and cutoff parameters are given in the previous subsection. Explicit expression 
for all these operators can be found in Ref. ||. Here we only note that: i) the tt- and p- 
meson exchange charge operators, the former of which gives by far the dominant contri- 
bution, are constructed using the PS (7r-like) and V (p-like) components projected out of 
the isospin-dependent spin-spin and tensor terms of the interaction ||, thus reducing their 
model dependence. The resulting two-body operators are denoted as PS and V, and are here 
obtained from the charge-independent part of the AV18. ii) In the pion (as well as vector 
meson) charge operators there are additional contributions due to the energy dependence 
of the pion propagator and direct coupling of the photon to the exchanged pion (p-meson). 
However, these operators give rise to non-local isovector contributions which are expected 
to provide only small corrections to the leading local terms. For example these operators 
would only contribute to the isovector combination of the 3 He and 3 H charge form factors, 
which is anyway a factor of three smaller than the isoscalar. Thus they are neglected in the 
present work. 



C. The three-body exchange current associated with S'-wave pion rescattering 

The isospin odd "large" component of the S'-wave pion-nucleon (ttN) scattering ampli 
tude at low energy and momentum transfer may be described by the effective interaction [|29| 



1 



AnrAW = -77^7 M T ' tp<t> X • (2-16) 



s 



Here <f> is the isovector pion field and /„• the pion decay constant (~ 93 MeV). This effective 
Lagrangian implies the "Weinberg- Tomozawa" relation for the isospin odd combination of 
the irN S -wave scattering lengths o,\ : 



A 3 
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m 



(ax - a 3 ) 



1 ( m m 



16tt V f„ 



(2.17) 



which agrees well with the experimental scattering length values. Combined with the pseu- 
dovector ttNN interaction 



NN 



firNN 



(2.18) 



where /, 



-y.y — 1, this interaction gives rise to the three-body interaction: 

v 2 



V s = — 



NN 



Am $1 \ 



Ti ■ Tj X T fc 

ijk cyclic 



& , ■ ki<T fc ■ k fc 
D~Dk 



l r 



o-j ■ ki x k fc + 



ki ■ [(Pi + Pi) - (Pi + Pj)} - k fc ■ [(p fc + p' k ) - ( Pj + p',)] • (2.19) 



The momentum vectors are defined so that kj denotes the fractional momentum transfer to 
nucleon i. The denominator factors D; are defined as 



Di = k + mi 



[2.20) 



The derivative couplings in the Lagrangians ( |2.16| ) and ( 2.19 ) lead to electromagnetic 
contact terms. These may be constructed by minimal substitution, and are found to have 
the expressions 



C 



TTTT-yNN 



(2.21) 



C^nn = -^^7V^(r x d>) z i> , (2.22) 

m n 

respectively. When complemented with the electromagnetic coupling of the pion, 

= -A^cf> x d»cj>) z , (2.23) 

these contact terms give rise to the following set of three-nucleon exchange current operators: 
(a) a contact current at the S'-wave rescattering vertex, (b) two contact currents at the two 
accompanying pseudovector irNN vertices and (c) two pion current terms. The explicit 
expressions for these are in the corresponding order: 

j ^' fc(q) = 8^7| fef) [Tfc x (Tj x Ti) + TiX (Tj x Tk ^ z 

pppj '-fa x (q - k, - kfc) - i(Pj + p ■) , (2.24) 
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[Ti X (-7- X T fc )] 



(tAcfu • let 



[<r,- • (ki -q)xk t 



l r 

+ 2 



4m /2 y 

k t ■ [(Pi + Pi) - (Pj + Pj)} - k fe ■ [(p fc + p k ) - (pj + p'j)} 



2m u + q • (pj + p'j) 



k) 



(2.25) 



i 1 ( Unn\" < ( \ i (<7i ■ ki)(cr fe • kjfc) 2k { - q 



4m /2 V m^ 



l r 



[ctj • (ki - q) x k fc ] + - [ki ■ [(pi + p 4 ) - ( Pi + Pj )] 



- k fc ■ [( Pfe + P ^) - (pj + p,)] - 2mw + q • (pj + p' ? 



+ i 



(2.26) 



In these exchange current operators the fractions of the total momentum transfer q imparted 
to the three nucleons are denoted k, respectively so that q = ki + k.2 + k3 . The denominator 
factors Di are defined in (|2.20|) , while the denominator factors D/ are defined as 



k,) 2 + ml 



(2.27) 



The combined three-nucleon exchange current operator j a + j b + j c satifies the continuity 
equation with the three-nucleon interaction Vs (|2.19|) , as may be verified by comparing the 
product q • j with the commutator of Vs and the single nucleon charge operator. These 
two-pion exchange three-nucleon currents will be labelled as tttts below. 

Note that the three-nucleon interaction (|2.19 ) is not contained in the Urbana IX three- 
nucleon interaction model [Q, the main part of which takes into account exchanges that 
involve excitation of intermediate A-isobar resonances, which are treated explicitly below. 
It should however be included in any complete three-nucleon interaction model, as it is 
implied by effective chiral Lagrangian models for the pion-nucleon system. It is included in 
three-nucleon interaction models that are based on pion exchange and rescattering described 
by current algebra or chiral Lagrangians (cf in the "d" term of the three-nucleon interaction 
in Ref. Bg). 



D. Elastic form factors of 3 H and 3 He 

In this section we present results for the magnetic moments, charge and magnetic form 
factors of 3 H and 3 He. The nuclear ground states are described by the PHH wave function 
obtained from the AV18 /UIX Hamiltonian model. 

A convenient espression to calculate the magnetic form factors of a J =1/2 nucleus, such 
as the A=3 systems under consideration here, is obtained by orienting the coordinate system 
so that the spin-quantization axis (the z-axis) lies along the momentum transfer q. It is 
then found that 

2m 1 

F M (q) = (*+|j x (gz)|tt_) , (2.28) 
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where \x is the nuclear magnetic moment in terms of /ijv (the nuclear magneton), \?+/_ are 
the normalized ground state wave functions with J z = ±1/2, respectively, and j x (qz) is the 
x-component of the current operator. Note that Fm{0) =1. The charge form factor is easily 
obtained from 

F c (q) = ^(^ + \p(qz)\^ + ) , (2.29) 

with F c (0) = 1. 

The matrix elements ( |2.28| ) and ( 2.29|) are evaluated with Monte Carlo methods. The 



wave function is written as a vector in the spin-isospin space of the three nucleons for any 
given spatial configuration R = (ri,r 2 ,r 3 ). For the given R, the state vectors j x (qz) | 
and p(qz) \ *ff + ) are calculated by performing exactly the spin-isospin algebra with the tech- 
niques described in Refs. HH- The spatial integrations are carried out by sampling the 
R-configurations according to the Metropolis et al. algorithm |5T]]. Typically, 400,000 con- 
figurations are enough, in the form factor calculations reported here, to achieve a relative 
error of few % at low and moderate values of momentum transfer q (q < 5 fm _1 ), increasing 
to ~30% at the highest g-values. 



1. The magnetic form factors 



The current operator includes, in addition to the one-body current in Eq. ( |2.13|) , the 
MI two-body currents obtained from the charge-independent part of the AV18 interaction 
(denoted as PS or 7r-like, V or p-like, SO, LL and S02), the MD pir'y and uotv) two-body 
currents, and finally the local terms of the three-body current associated with the S"-wave 
two-pion exchange three-nucleon interaction (|2.19| ) described in the previous section. 

Because of destructive interference in the matrix element for the magnetic dipole 
transition between the S- and D-state components of the wave function, the one-body 
predictions for the 3 H and 3 He magnetic form factors (MFF) have distinct minima at 
around ~3.5 fm _1 and ~2.5 fm^ 1 , respectively, in disagreement with the experimental 
data |3^j3^j3^j3^j3^j3T|j3^j3^j40|j4l| , as shown in Figs. [I] and |2|. The situation is closely 
related to that of the backward cross section for electrodisintegration of the deuteron, which 
is in fact dominated by two-body current contributions for values of momentum transfer 
above ~2.5 fm" 1 |^2| . 



Inclusion of the contributions from the two- and mis three-body currents shifts the zeros 
in the calculated MFF to higher g-values. While the experimental 3 H MFF is in good 
agreement with theory over a wide range of momentum transfers, there is a significant 
discrepancy between the measured and calculated values of the 3 He MFF in the region of 
diffraction minimum. This discrepancy persists even when different parametrizations of the 
nucleon electromagnetic form factors are used. This is evident from Figs. [I] and |2] where the 
total results obtained with the Gari-Krumpelmann (GK) parametrization |Q of the nucleon 
electromagnetic form factors are shown. 

It is useful to define the quantities 



F S ^(q) 



/i( 3 He)F M (g; 3 He) ± /i( 3 H)F M (g; 3 H) 



^2.30) 
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If the 3 H and 3 He ground states were pure T—l/2 states, then and linear combina- 
tions of the three-nucleon MFF would be only influenced by, respectively, the isoscalar (S) 
and isovector (V) parts of the current operator. However, small isospin admixtures with 
T >l/2, induced by electromagnetic, CSB and CIB terms present in the AV18 interaction, 
are included in the present wave functions. As a consequence, purely isoscalar (isovector) 
current operators give small, otherwise vanishing, contributions to the F^ (F^ d ) MFF. 

The contributions of the individual components of the two- and three-nucleon (nns term) 
currents to the F^ and F^ combinations are shown in Figs. |3] and |j. In the diffraction 
region the PS (7r-like) isovector current gives the dominant contribution to F^, while the 
contributions from remaining currents are significantly smaller, about one order of magnitude 
or more. The three-nucleon current (titts) associated with the S"-wave irN coupling is found 
to give a very small correction. 

Among the two-body contributions to F^, the most important is that due to the cur- 
rents from the spin-orbit interactions (SO), the remaining operators producing a very small 
correction. Note that the isovector PS and V currents contribute to F^ because of the small 
isospin-symmetry breaking components present in the 3 H and 3 He wave functions induced 
by the AV18 model, as mentioned earlier. 

Finally, the cumulative contributions to the F^ and F^ combinations are compared with 
the experimental data [41] in Figs. [5] and |6], respectively. The zero in the calculated F^ is 
found to occur at lower q- value than experimentally observed. As shown in the next section, 
this discrepancy between theory and experiment remains unresolved even when A-isobar 
degrees of freedom are included in both the nuclear wave functions and currents. We will 
return to this point in the conclusions. Predictions for the magnetic moments are given in 
Tables |I| and III, while those for the magnetic radii are listed in Table IV. These results 
are discussed in Sec. III-D. 



2. The charge form factors 

The charge operator includes, in addition to the one-body term of Eq. (|2.12|) , the PS 
or 7r-like, V or p-like, u, pttj and unj two-body operators, discussed previously. The cal- 
culated 3 H and 3 He charge form factors (CFF) are compared with the experimental data 
^ , |33|j3^ , |35lJ3^ , |37|j38] , |39|^ , |4ll| in Figs. [7| and |8|. There is excellent agreement between theory 



and experiment, as is clear from these figures. The important role of the two-body contri- 
butions above 3 fm -1 is also evident. The remarkable success of the present picture based 
on non-relativistic wave functions and a charge operator including the leading relativistic 
corrections should be stressed. It suggests, in particular, that the present model for the two- 
body charge operator is better than one a priori should expect. These operators, such as 
the PS charge operator, fall into the class of relativistic corrections. Thus, evaluating their 
matrix elements with non-relativistic wave functions represents only the first approximation 
to a systematic reduction. A consistent treatment of these relativistic effects would require, 



for example, inclusion of the boost corrections on the nuclear wave functions p5| , p6| , ^3 
Yet, the excellent agreement between the calculated and measured CFF suggests that these 
corrections may be neglegible in the g-range explored so far. 

For completeness, we show in Figs. |9] and |10| the contributions from the individual com- 
ponents of the charge operator to the linear combinations 
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?s,v, 



c (q) = -[2F c (q; i Re)±F c (q;"R) 



(2.31) 



Note that again, because of isospin-symmetry breaking components present in the 3 He and 
3 H wave functions, the purely isovector (isoscalar) unr-f (pn'y) charge operator gives a small, 
otherwise vanishing, correction to the Fq (Fq) CFF. 

Finally, values for the charge radii of 3 H and 3 He are listed in Table |V[ The results 
including the contributions associated with the two-body charge operators are found to be 
in good agreement with experimental data. 



III. BEYOND NUCLEONS ONLY 

The simplest picture views the nucleus as being made up of nucleons, and assumes that 
all other sub-nucleonic degrees of freedom may be eliminated in favor of effective many-body 
operators acting on the nucleons' coordinates. The validity of such a description is based on 
the success it has achieved in the quantitative prediction of many nuclear observables [[!]. 
However, it is interesting to consider corrections to this picture by including the degrees of 
freedom associated with nuclear resonances as additional constituents of the nucleus. When 
treating phenomena which do not involve explicitely meson production, it is reasonable to 
expect that the lowest excitation of the nucleon, the A-isobar, plays a leading role. 

In such an approach, the A=3 nuclear wave function is written as 

\$j n+a = ^(NNN) + ¥ 1] (NNA) + ¥ 2) (NAA) + ^ (3) (AAA) , (3.1) 

where \1/ is that part of the total wave function consisting only of nucleons; the term 

is the component in which a single nucleon has been converted into a A-isobar, and so on. 

The nuclear two-body interaction is taken as 



Vij 



£ £ B\b\) , (3.2) 

B^B^NA B'B'=N,A 



where transition interactions such as Vij(NN — > NA), Vij(NN — > AA), etc. are responsible 
for generating A-isobar admixtures in the wave function. The long-range part of is 
due to pion-exchange, while its short- and intermediate-range parts, influenced by more 
complex dynamics, is constrained by fitting NN scattering data at lab energy < 400 MeV 
and deuteron properties 

Once the NN, NA and AA interactions have been determined, the problem is reduced 
to solving the N-A coupled-channel Schrodinger equation. In principle, for the A=3 sys- 
tems Faddeev and hyperspherical-harmonics techniques can be used (and, indeed, Faddeev 
methods have been used in the past |45| , ^6| ) to this end, although the large number of 
N-A channels involved makes the practical implementation of these methods difficult. A 
somewhat simpler approach consists of a generalization of the correlation operator tech- 
nique |47| , which has proven very useful in the variational theory of light nuclei, particularly 
in the context of variational Monte Carlo calculations 0,HH[ . In such an approach, known as 



the transition-correlation-operator (TCO) method f4"9 ], the nuclear wave function is written 

as 
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N+A 



sl[(i + u, 



TR 



i<j 



(3.3) 



where \I/ is the purely nucleonic component, S is the symmetrizer and the transition operators 
U$ R convert NN pairs into AA and AA pairs. In the present study the ^ is taken from 
PHH solutions of the AV18/UIX Hamiltonian with nucleons only interactions [j|9], while the 
Ufj R is obtained from two-body bound and low-energy scattering state solutions of the full 
iV-A coupled-channel problem. This aspect of the present calculations is reviewed briefly in 
the next subsection. 



A. Wave functions with A-admixtures 



The transition correlation operator (TCO) method [415 consists in approximating the 
*$?n+a as in Eq. ( |3.3|) , with the transition operators Uj R defined as 



U 



TR 

ij 



(3.4) 



U NA 



U, 



ij 



rll. 



J J 



u 



(3.5) 
(3.6) 



Here, S« and Tj are spin- and isospin-transition operators which convert nucleon i into a 
A-isobar; S-j and S^ 11 are tensor operators in which, respectively, the Pauli spin operators 
of either particle i or j, and both particles % and j are replaced by corresponding spin- 
transition operators. The Vj R vanishes in the limit of large interparticle separations, since 
no A-components can exist asymptotically. 

The transition operator Uj R and nucleonic wave function \l> in Eq. ( |3.3|) could be de- 



termined variationally by using an interaction of the form given in Eq. ( |3.2| ), that contains 
both A^ and A degrees of freedom, such as the Argonne v 2 sq (AV28Q) model f?|j50|, and by 
minimizing the ground-state energy of each given nucleus. Instead, we use transition cor- 
relation functions u aTlI (r), etc. (shown in Fig. |ll|) that approximately reproduce two-body 
bound- and low-energy scattering-state wave functions for the AV28Q model, and take the 
PHH wave function obtained in Sec. II-A as the \l/ in Eq. (|3.3|). The validity of such an 



approximation has been discussed at length in the original reference |^9| . Here, we only note 
that i) since the correlation functions M <jr// (r), etc. are short-ranged, they are expected to 
have a rather weak dependence on A; ii) it is important the \I> used in Eq. ( |3.3|) , obtained 



from a Vij(NN — > NN) interaction phase-equivalent to the full Vij of Eq. ( |3.2| ), be propor- 
tional to that projected out from the ^tv+a wave function for the interaction. This has 
been explicitly verified by direct calculation in the two-body problem HDI . 

In the TCO scheme, the perturbation theory description of A-admixtures is equivalent 
to the replacements: 



U. 



VA,PT 



Vij(NN -> NA) 
m — ira 
^aa.pt = VyjNN -> A A) 
13 2(m — m/±) 



(3.7) 
(3.8) 
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where the kinetic energy contributions in the denominators of Eqs. (|3.7|) and Q3.8Q have 
been neglected (static A approximation). The transition interactions Vij(NN — > NA) and 
Vij(NN — > A A) have the same operator structure as U?f A and U AA of Eqs. ( |3.5| ) and Q3.6|), 
but with the u <TTa {r) and u tTa (r) functions replaced by, respectively, 

v™( r ) = ^^^C(x), (3.9) 
«*"(r) = (l + * + A) — C\x) . (3.10) 

47T 3 V X X z J X 

Here a = II, III, x = m/, (//) a = UnnUna, Una/kMA for a = II, III, respectively, and 
the cutoff function C(x) = 1 - e~ A:c2 , A = 4.09 in the AV28Q model Q. Note that in 
Fig. O u qll,FT (r) = v gII (r)/(m — m&) and u qIll,FT (r) = v qlu (r)/[2(m — m^)], with q = ar, 



tr. 

The perturbative treatment has been often (in fact, almost exclusively) used in the 
literature to estimate the effect of A degrees of freedom on electroweak observables. However, 
it may lead to a substantial overprediction of their importance [0,0 , since it produces NA 



and A A wave functions which are too large at short distance, see Fig. 11 



B. iVA-transition and A currents 

The nuclear electromagnetic current is expanded into a sum of many-body terms that 
operate on the nucleon and A-isobar degrees of freedom. The nucleonic component of this 
current operator has been discussed in the previous section. Here, we only discuss its A 
components. 



1. N A-transition and A one-body currents 
The one-body current is written as 

i (1) (q)= E j*(qjfl-s'), (3-ii) 

b,b'=n,a 

where j«(q; — > N) is the nucleonic current component given in Eq. ( J2.13 ) and 



j,(q; N - A) = -^G 7M (gJ)^q x S,T 2ji , (3.12) 
j,(q; A -+ A) = -^-G TAA (gJ)e iq ' ri q x E<(1 + Q z>i ) . (3.13) 

Here S (0) is the Pauli operator for the A spin 3/2 (isospin 3/2), and the expression for 
ji(q; A — > N) is obtained from that for jj(q; N —>■ A) by replacing the transition spin and 
isospin operators by their hermitian conjugates. The N A-transition and A electromagnetic 
form factors, respectively G 7 tva and G 7 aa, are parametrized as 
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G 7 aa(^) 



( 



fJ'-yAA 



(l + ^/A 2 AA 



(3.14) 
(3.15) 



Here the iVA-transition magnetic moment ^na is taken equal to 3 /ijv, as obtained from 
an analysis of 7 AT data in the A-resonance region [51|; this analysis also gives AjvA,i — 
0.84 GeV and A^a.2 = 1-2 GeV. The value used for the A magnetic moment /i 7 AA is 
4.35 /ijv by averaging results of a soft-photon analysis of pion-proton bremsstrahlung data 
near the A ++ resonance | 52|| , and Aaa = 0.84 GeV as in the dipole parametrization of 
the nucleon form factor. In principle, N to A excitation can also occur via an electric 
quadrupole transition. Its contribution, however, has been ignored, since the associated 
pion photoproduction amplitude is found to be experimentally small at resonance [53]. Also 
neglected is the A convection current. 



2. N A-transition two-body currents 



The two-body term is written as 



jl?(q) 



b' b'.=n,a j^'( q ' BiB i 



(3.16) 



where the prime over the summation symbols indicates that terms involving more than a 
single A have been neglected in the present study. The NN —>■ NN two-body terms have 
already been discussed. The two-body terms involving at most a single A are illustrated in 
Fig. [12], and are explicitly given by 



jii(q; NN - NA) = (n x T 



3)z 



+ ( *i'X*i(tr i . V i )(S i - V^hir, 



101 



(3.17) 



where = 
respectively, 



Tj, = (rj + Tj)/2, and the functions h(r) and h(r) are defined as, 



h(r) 
h(r) 



Arc 



firNNfiT 



NA 



An 



At / ' 2 dze-' lz * r e~ rL{z) 
mt 



(3.18) 
(3.19) 



with x = m^r and L(z) = [m^ + q 2 (l/A — z 2 )] 1 ^ 2 . Terms explicitly proportional to q in 
Eq. (|3.17 ) have been dropped, since in applications only the transverse components of j(q) 
occur. The three terms in Eq. (|3.17|) are associated with diagrams (a), (b) and (c) in Fig. [l^, 
respectively, and can be obtained from the well known expression of the two-body nucleonic 
currents due to pion-exchange by replacing a-j and Tj with Sj and Tj, respectively. 
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To account for the hadron extended structure, form factors must be introduced at the 
ttNN and ttNA vertices. In the case of Vij(NN — > iVA) interaction, an r-space gaussian 
cutoff has been used. However, for the j(NN — > NA) above it has been found convenient to 
introduce monopole form factors with A=900 MeV in its p-space expression. This value for 
A is consistent with that obtained from the tensor component of Vij(NN — > NA). Finally, 
the expression in Eq. (|3.17 ) is multiplied by the isovector form factor G\(qfy. 



C. Calculation 

Calculation of the magnetic form factors requires evaluation of the transition matrix 
element in Eq. ( 2.28Q , where the wave functions and currents include both nucleonic and 



A-isobar degrees of freedom. To evaluate such a matrix element, it is convenient to expand 
the wave function ^n+a,j z as 

*n + a,j z = *j z + E Vf**J. + • • • > (3.20) 



and write the numerator of Eq. (|2.28 ), in a schematic notation, as 



(*N+Aj\j\*N + A,i) = (tf/|j(Wonly)|*i) + (^ f \j(A)\^ i ) , (3.21) 

where j(N only) denotes all one- and two-body contributions to j(q) which only involve 
nucleon degrees of freedom, i.e., j(N only) = j<$(N -> N) + j®(NN -> NN). The 
operator j(A) includes terms involving the A-isobar degrees of freedom, associated with the 
explicit A currents j {1) (N -> A), j (1) (A -> iV), j {1) (A -> A), (NN -> NA), etc., and 
with the transition operators lfE R . The operator j(A) is illustrated diagrammatically in 



Figs. |13| and The terms (a)-(g) in Fig. ^ and (a)-(f) in Fig. |14| are two-body current 



operators. The terms (g)-(l) in Fig. are three-body current operators, while the terms (h)- 
(j) in Fig. 13 are to be interpreted as renormalization corrections to the "nucleonic" matrix 
elements | j(A^only) | \I/A, due to the presence of A-admixtures in the wave functions. 

There are, however, additional, connected three-body terms in j(A) that are neglected 
in the present work. A number of these are illustrated in Fig. [15]. Their contribution is 



expected to be significantly smaller than that from the terms in Figs. 13 and O involving 

,\ — ; — 

transition correlations between two particles only, of the type U^ B U BB , but comparable 



to that from the three-body terms in Fig. [14| having U^ B UK 3 . These have been found to 



'ij u jk 

be very small. 

The terms in Fig. [13| are expanded as operators acting on the nucleons' coordinates. For 
example, the terms (a) and (e) in Fig. [T3] have the structure, respectively, 



(*)=jl 1 \A-+N)U$ N , (3-22) 
(e) = ^Sf (A^A)C/^, (3.23) 



which can be reduced to operators involving only Pauli spin and isospin matrices by using 
the identities 
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S t -AS-B = ^A-B-^cr-(AxB) , (3.24) 

S t -AE-BS-C = -iA-(BxC)--er.AB-C 

-~A BC <r + ~A (B <r)C , (3.25) 
3 3 

where A, B and C are vector operators that commute with cr, but not necessarily among 
themselves. 



While the terms in Fig. [14] could have been reduced in precisely the same way, the 
resulting expressions in terms of cr and r Pauli matrices become too cumbersome. Thus, 
for these it was found to be more convenient to retain the explicit representation of S (S^) 
as a 4 x 2 (2 x 4) matrix 

/ -e_ \ 

2 eo -^?e_ 



V -e+ J 



where e± = =f(x ± iy)/y/2, eo = z, and e* = (— )^e_ M and derive the result of terms such 

as ( a ) + (c) + (e)=[7^ At j\f{NN -> NA) on the state by first operating with and 

then with U NA \ The Monte Carlo evaluation of the matrix element is then performed with 
methods similar to those sketched in Sec. II-D. 

The normalization of the wave function is given by 

(*n + a,j z | *n + a,j,) = 1 1 + E[2o£"V^ + u^ ] u^\ i * Jb > 

i<j 

+ ( three — body terms) (3.26) 

and the three-body terms have been neglected consistently with the approximation intro- 
duced in Eq. ( |3.21|) , as discussed above. 

Perturbation theory (PT) estimates of the importance of A-isobar degrees of freedom 
in photo- and electro-nuclear observables typically include only the contribution from single 
iV # A transitions (namely diagrams (a) and (b) in Fig. [13]) and ignore the change in the 
wave function normalization. In the TCO scheme, these PT estimates are obtained by using 
jjBB ,pt transition correlation defined in Eqs. ( ]3.7|) and (|3.8j ) ||49|| . In particular, the PT 
expressions for the three-body terms in Fig. |T4], diagrams (g)-(h)along with those in which 
the first and third legs are exchanged, can easily be shown to satisfy current conservation 
with the Fujita-Miyazawa two-pion exchange three- nucleon interaction (TNI) given by 

V? k M (NNN -> NNN) = v ik (AN -> NN) ViANN -> NA) + h.c. , (3.27) 

J m — mA 



where the transition potentials are defined in Eqs. (|3.9| ) and ( |3.10| ) (here with the cutoff 



function C(x) set to one). Current models of TNI M include the "long-range" 27r-exchange 
component above. Indeed, the need of including the associated three-body currents provided 
one of the motivations for undertaking the present study. 
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D. The magnetic moments and form factors 



The 3 H and 3 He magnetic form factors obtained by including nucleon and A-isobar 



degrees of freedom in the nuclear wave functions and currents are shown in Figs. 16 and 



17; while individual contributions to the combinations Fj^ and FY, are displayed in Figs. |18 



M aLlu - ± M 

and [19]. Finally, individual and cumulative contributions to the magnetic moments and 
cumulative contributions to the magnetic radii of the trinucleons are listed in Tables |T|, [1 1 1 
and [IV], respectively. Note that in Figs. [16] and [H] and Table [TT] the contributions labelled 



1-A and 2-A are associated with the diagrams in Figs. [13| and [LJ, respectively. Also note 
that the individual nucleonic and A-isobar contributions in Figs. |TB] and [TP| and Table |TJ are 
normalized as, in a schematic notation, 



However, the cumulative contributions in Figs. [16] and [17] and Table |T| and |TV] are normal- 
ized as 

[TOT— N] = <*'^y*> , (3.29, 
when "nucleons only" terms are retained, and as 



(^N+A 


j(N + A) 


^n+a) 




Vn+a) 



[TOT- (N + A)] = "^'^ — 71 iV ^ / , (3.30) 



when, in addition, the A terms are included. 

The contributions associated with A-components are found to be small in contrast to 
earlier studies . In particular, we find that the sign of the 2-A contribution in Fig. |TU 



is 



opposite to that reported in Ref. The origin of this difference is unclear at this point. 



However, we do find that the sign of the 2-A contribution (see Fig. |T^) is the same as that 
of the nucleonic PS (7r-like) contribution, as one would expect. 

The predicted magnetic moments of the trinucleons are within less than 1% of the exper- 
imental values. The predominatly isovector A-isobar contributions lead to an increase (in 
magnitude) of the 3 H and 3 He magnetic moments calculated with nucleons only degrees of 
freedom of, respectively, 1.1% and 1.7% in relative terms. We note that perturbation theory 
estimates of the A-isobar contributions are found to be significantly larger than obtained 
here [[|. 

The predicted magnetic radii of 3 H and 3 He are, respectively, 2% and 3% smaller than 
the experimental values, but still within experimental errors. Inclusion of the contributions 
due to two- and three-body exchange currents leads to a decrease of the 3 H and 3 He magnetic 
radii of, respectively, 5% and 6%. 

While the agreement between theory and experiment is satisfactory for the magnetic 
moments, magnetic radii and low q form factors, the calculated form factors, particularly 
that of 3 He, remain at variance with the experiment in the diffraction region. The role played 
by A-isobar degrees of freedom is found to be marginal over the whole g-range considered 
here. 
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IV. CONCLUSIONS 



The present results for the electromagnetic form factors of the trinucleons may be sum- 
marized as follows: i) the trinucleon charge form factors agree well with the experimental 
values when calculated with wave functions obtained from a Hamiltonian consisting of the 
Argonne v\$ two-nucleon and the Urbana IX three-nucleon interactions; ii) agreement with 
the experimental charge form factors requires that the two-nucleon exchange charge opera- 
tors are taken into account; iii) the calculated magnetic form factor of 3 H agrees well with 
experiment, whereas that of 3 He agrees well with the experimental values only for momen- 
tum transfer values below the first zero in the form factor; iv) the two-nucleon exchange 
current contributions are essential for achieving agreement with experiment whereas three- 
nucleon exchange current operators and the A-isobar configurations have only very small 
effects on the calculated magnetic form factors. 

The result for the static observables are that the calculated value for the isovector com- 
bination of the trinucleon magnetic moments agrees completely with the experimental value 
(Table |ITJ). The isoscalar combination of the trinucleon magnetic moments exceeds the 
experimental value by about 5%, but this small disagreement does not prevent a good re- 
production of the isoscalar combination of the experimental magnetic form factors. As the 
calculated magnetic moments of 3 H and 3 He differ by less than 0.015 n.m. from their exper- 
imental values, the results appear to be very satisfactory. The calculated charge radii are 
smaller by only 2% than the experimental values. The calculated magnetic radii are smaller 
than 3% than the experimental values (Tables |V| and 0). To obtain these quite satisfactory 
calculated values for the charge and magnetic radii the exchange current contributions have 
to be taken into account. 

We note finally that the three-nucleon exchange current operator (|2.24j) - (|2.26|) , which 
was constructed to satisfy the continuity equation with the three-nucleon interaction ( [2.19P , 
was found to give only very small contributions to the magnetic form factors of the trin- 
ucleons. It is worth noting that the corresponding component of the "Tucson-Melbourne" 



type three-nucleon interaction in |30[ is roughly an order of magnitude stronger than the 
three-nucleon interaction ( |2.19p , and would therefore imply correspondingly much larger 
three-nucleon exchange current contributions. Part of this difference is the inclusion of 
A-isobar intermediate state effects in the "d" term of that three-nucleon interaction. In 
the present work A-isobar configurations are treated explicitly, and should therefore not be 
included in irreducible three-nucleon exchange current operators. 
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TABLES 



TABLE I. Binding energies corresponding to the AV14 and AV18/UIX Hamiltonian models. 
The AV14 results obtained with the PHH expansion are compared with those calculated by solving 
the Faddeev equations in configuration (F/R) and in momentum (F/P) space. The statistical error 
associated with the GFMC calculations are shown in parenthesis. 



Model 


Method 


B( 3 H) (MeV) 


B( 3 He) (MeV) 




PHH 


7.683 




7.032 


AV14 


F/R 


7.670 




7.014 




F/P 


7.680 






AV18/UIX 


PHH 


8.49 




7.75 




GFMC 


8.47(1) 




7.71(1) 


expt. 




8.48 




7.72 


TABLE II. Individual contributions from the different components of the nuclear electromag- 


netic current operator to the 3 H and 3 


He magnetic moments and their /xg and fly 


combinations, in 


nuclear magnetons (n.nL 


). Note that, because of isospin-symmetry breaking components present in 


the PHH 3 H and 3 He wave functions, 


purely isoscalar (isovector) currents give non vanishing con- 


tributions to the fly (l^S, 


) combination. The contributions to ^5 


due to the tttts 


and 2-A currents 


and those to nv due to the S02+LL currents are very small and are not listed. 








M 3 He) 


MS 


Hv 


1-N 


2.571 


-1.757 


0.407 


2.164 


PS 


0.274 


-0.269 


0.002 


0.271 


V 


0.046 


-0.044 


0.001 


0.045 


SO 


0.057 


0.010 


0.033 


0.023 


S02+LL 


-0.005 


-0.006 


-0.005 






0.016 


-0.009 


0.003 


0.012 


TTTTs 


0.002 


-0.002 




0.002 


1-A 


0.084 


-0.064 


0.010 


0.074 


2-A 


0.024 


-0.024 




0.024 


TABLE III. Cumulative and normalized contributions to the 3 H and 3 He magnetic moments 


and their fis and \xy combinations, in 


nuclear magnetons (n.m.). 


compared with the experimental 


data. 












M 3 H) 


^( 3 He) 




Hv 


1-N 


2.571 


-1.757 


0.407 


2.164 


TOT-N 


2.961 


-2.077 


0.442 


2.519 


TOT-N+1-A 


2.971 


-2.089 


0.441 


2.530 


TOT-(N+A) 


2.994 


-2.112 


0.441 


2.553 


EXP 


2.979 


-2.127 


0.426 


2.553 
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TABLE IV. Cumulative and normalized contributions to the 3 H and 3 He r.m.s. magnetic radii, 
in fm, compared with the experimental data. 





3 H 


3 He 


1-N 


1.895 


2.040 


TOT-N 


1.810 


1.925 


TOT-N+1-A 


1.804 


1.916 


TOT-(N+A) 


1.800 


i .yuy 


EXP 


1.840±0.181 


l.yboitu.loo 


TABLE V. Cumulative and normalized contributions to the 3 H and 3 He r.m.s. 


charge radii, in 


fm, compared with 


the experimental data. 






3 H 


3 He 


1-N 


1.711 


1.919 


TOT 


1.725 


1.928 


EXP 


1.755±0.086 


1.959±0.030 
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FIGURES 



FIG. 1. The magnetic form factors of 3 H, obtained with single- nucleon currents (1-N), and 
with inclusion of two-nucleon current [(l+2)-N] and irirs three-nucleon [TOT-N(D)] current con- 
tributions, are compared with data (shaded area) from Amroun et al. [El]]. Theoretical results 
correspond to the AV18/UIX PHH wave functions, and employ the dipole parametrization (in- 
cluding the Galster factor for Ge[<&)) for the nucleon electromagnetic form factors. Note that the 
Sachs form factor G^(g^) is used in the model-independent isovector two-body currents obtained 
from the charge-independent part of the AV18 interaction. Also shown are the total results cor- 
responding to the Gari-Krumpelmann parametrization 43] of the nucleon electromagnetic form 
factor [TOT-N(GK)]. 



FIG. 2. Same as in Fig. @, but for 3 He. 

FIG. 3. Individual contributions to the F^icta) combination, Eq. (|2.30| ), of the 3 H and 3 He 
magnetic form factors, obtained with the dipole parametrization of the nucleon electromagnetic 
form factors. The sign of each contribution is given in parenthesis. Note that, because of 
isospin-symmetry breaking components present in the 3 H and 3 He wave functions, the purely 
isovector PS, V and irirs currents give non vanishing contributions to the i^(^) combination. 
However as the 7r7r,s contribution is very small, is not shown. 



FIG. 4. Individual contributions to the F^(q ll ) combination, Eq. ( 2.30| ), of the 3 H and 3 He 



magnetic form factors, obtained with the dipole parametrization of the nucleon electromagnetic 
form factors. The sign of each contribution is given in parenthesis. Note that, because of 
isospin-symmetry breaking components present in the 3 H and 3 He wave functions, the purely 
isoscalar pir-y current gives non vanishing contributions to the F^ I (q fl ) combination. However, 
being very small, it is not shown. 

FIG. 5. The F^, I (q fl ) combinations of the 3 H and 3 He magnetic form factors, obtained with 
single-nucleon currents (1-N), and with inclusion of two-nucleon current [(l+2)-N] and irirs 
three-nucleon current (TOT-N) contributions, are compared with data (shaded area) from Amroun 
et al. [El]]. The dipole parametrization is used for the nucleon electromagnetic form factors. 



FIG. 6. Same as in Fig. |5], but for the F^(q^) combination of the 3 H and 3 He magnetic form 
factors. 

FIG. 7. The charge form factors of 3 H, obtained with a single-nucleon charge operator (1-N) 
and with inclusion of two-nucleon charge operator contributions (TOT-N), are compared with data 



(shaded area) from Amroun et al. 41]. Note that the 1-N results also include the Darwin- Foldy 
and spin-orbit corrections. Theoretical results correspond to the AV18/UIX PHH wave functions, 
and employ the dipole parametrization of the nucleon electromagnetic form factors. 
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FIG. 8. Same as in Fig. 0, but for 3 He. 



FIG. 9. Individual contributions to the F^q^) combination, Eq. ( 2.31 ), of the 3 H and 3 He 



3xj o^^l 3t 



charge form factors, obtained with the dipole parametrization of the nucleon electromagnetic 
form factors. The sign of each combination is given in parenthesis. Note that, because of 
isospin-symmetry breaking components present in the 3 H and 3 He wave functions, the purely 
isovector ujir^f charge operator gives a non vanishing contribution to the F^q^) combination. 



FIG. 10. Individual contributions to the F^{q fl ) combination, Eq. ( 2.31| ), of the 3 H and 
3 He charge form factors, obtained with the dipole parametrization of the nucleon electromag- 
netic form factors. The sign of each combination is given in parenthesis. Note that, because 
of isospin-symmetry breaking components present in the 3 H and 3 He wave functions, the purely 
isoscalar ptrj charge operator gives a non vanishing contribution to the F^ (q^) combination. 

FIG. 11. Transition correlation functions u aTlI (r), u trII (r), etc. obtained for the AV28Q 
model [5C], and perturbation theory equivalents u aTlI,PT (r), u tTlI,PT (r), etc. 



FIG. 12. iVA-transition two-body currents due to pion exchange. 

FIG. 13. Diagrammatic representation of operators included in j(A) due to one-body currents 

jW(N -» A), jW(A -» N) and jW(A -> A), and transition correlations U NA , U AN , U AA , 

and corresponding hermitian conjugates. Wavy, thin, thick, dashed and cross-dashed lines denote 

' 't 

photons, nucleons, A-isobars and transition correlations U BB and U BB , respectively. 

FIG. 14. Diagrammatic representation of operators included in j(A) due to two-body currents 
j( 2 )(NN -> NA), j@)(NN -> AN), etc., and transition correlations U NA , U AN , and corre- 
sponding hermitian conjugates. Wavy, thin, thick, dashed and cross-dashed lines denote photons, 

nucleons, A-isobars and transition correlations U BB and U BB , respectively. 

FIG. 15. Diagrams associated with connected three-body terms, which are neglected in the 
present work. Wavy, thin, thick, dashed, cross-dashed and dotted lines denote photons, nucleons, 

A-isobars, transition correlations U BB and U BB , and the two-body curr ent jW(NN -» NN), 
respectively. 

FIG. 16. The magnetic form factors of 3 H, obtained with single-nucleon currents (1-N), and 
with inclusion of two- and three-nucleon current (TOT-N) and A [TOT-(N+A)] contributions. 



FIG. 17. Same as in Fig. 16, but for 3 He. 
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FIG. 18. The single-nucleon contribution to the F^^q^) combination of the 3 H and 3 He mag- 
netic form factors is compared with the 1-A and 2-A contributions, associated respectively with 
diagrams of Fig. [l^ and [l4|. 



FIG. 19. The single-nucleon and leading PS two-nucleon contributions to the F^ I (q fl ) combina- 
tion of the 3 H and 3 He magnetic form factors are compared with the 1-A and 2-A contributions, 
associated respectively with diagrams of Fig. O and [14. 
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